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1  Nomenclature 

A  =  b- 

s 

b/2 

aspect  ratio  of  the  wing; 

semispan  of  the  wing; 

Ci 

mean  chord  of  the  i  -th  part  of  the  wing; 

c 

C1 

dimensionless  chord  of  rectangular  outer  wing  part 
(referenced  to  the  tip  chord  of  the  main  part  of  the  wing); 

Csec 

section  chord; 

CA 

mean  aerodynamic  chord; 

Cm 

mean  chord; 

Ct 

tip  chord; 

Cr 

root  chord; 

c 

K 

K 

C*o 

c/ 

c 

ma 

l/dmax 

V  )  r 

b/2 

r  =  2  \r(z)dz 

0 

taper  ratio  of  the  main  part  of  the  wing  (referenced  to  the 
root  chord); 

slope  of  the  lift  curve; 

drag  due  to  lift  factor; 
leading-edge  sweepback  angle,  deg.; 

zero  lift  drag  coefficient; 

flat -plate  skin  friction  coefficient; 

derivative  of  a  pitching-moment  coefficient; 

maximum  aerodynamic  efficiency; 

dimensionless  velocity  circulation; 

net  velocity  circulation; 

Z 

spanwise  coordinate; 

Z  b/2 
t 

dimensionless  spanwise  coordinate; 

thickness  of  wing; 
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V 

M 


Vc 

Re,  =  -U 

v 

H 

v 

ma 


x c 


o 

N 

l 

R 


m 

n 

N 

Cp 

P 

P 

L 


m0 

mst 

Mowing 


dimensionless  airfoil  thickness; 

reference  wing  area,  sum  of  the  outer  (extendable)  part  area 
and  main  part  area. 

airspeed;  volume  of  structure  material; 

Mach  number; 

Reynolds  number  based  on  the  mean  chord  of  each  part  of 
the  wing; 
flight  altitude; 

coefficient  of  kinematic  viscosity,  m  /s; 

derivative  in  angle  atack  of  pitching  moment  coefficient; 

relative  mean  aerodynamic  center  of  wing. 

stress; 

internal  force; 
length; 

equivalent  force  flow; 

total  number  of  vortexes  on  the  wing; 

element  quantity; 

number  of  strips  on  a  semispan  of  the  wing; 
dimensionless  pressure  coefficient, 
density  of  structural  material; 
characteristic  (reference)  load; 
characteristic  (reference)  dimension  (length); 
take-off  mass; 
mass  of  structure; 
mass  of  wing  structure; 
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G 

Ck 
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New  criteria  and  notations  used  on  this  report: 

Load-Carrying  Factor  (LCF); 

dimensionless  coefficient  of  structure  load-carrying  perfection; 
dimensionless  full  mass  coefficient; 


Abbreviatures: 


FEM 

Finite  Element  Method; 

LCF 

Load-Carrying  Factor; 

TW 

Telescope  Wing; 

MW 

Morphing  Wing; 

MDV 

Method  of  Descrete  Vortexes; 

i 

eqv 

st 

T 


Indices: 

element  number, 
equivalent; 
structure; 
theoretical; 
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2  Summary 

The  present  report  contains  a  description  of  the  research  results  about  weight 
and  aerodynamic  efficiency  of  telescope  wings. 

Telescope  wings  (TW)  are  sufficiently  novel.  There  are  no  available  statistical 
data  for  such  wings.  Besides  the  efficiency  of  TW  should  obviously  depend  upon  a 
number  of  design  variables:  ratios  of  areas  and  lengths,  arrangement  of  stationary 
and  movable  parts  and  etc.  Known  weight  formulas  are  not  able  to  account  for  these 
features. 

That  is  why  in  the  present  work  a  vital  new  approach  for  calculation  of 
unconventional  aircraft  structures  is  proposed.  The  main  idea  of  this  approach,  when 
applied  to  TW,  is  the  following. 

1.  Into  the  geometrical  limits  of  the  TW,  including  internal  hollows,  a  three- 
dimensional  elastic  body  is  inscribed.  This  body  consists  of  some  hypothetical 
material  with  non-uniform  density.  We  consider  that  local  strength  of  this  material 
and  its  local  coefficient  of  elasticity  depend  linearly  upon  its  local  density. 

2.  Optimization  for  the  density  distribution  in  this  material  is  performed 
according  to  a  special  algorithm.  The  objective  function  for  this  optimization  if  the 
load-carrying  factor  (LCF)  G,  which  expresses  both  the  internal  loads  within  the 
structure  and  their  action  path  lengths.  In  case  of  3-D  body  it  is  calculated  as 

G  =  \oeqvdV 

V 

The  result  of  this  optimization  is  a  body  with  non-uniform  density  distribution 
(a  structure),  which  has  a  minimal  value  of  G  for  the  given  external  loading, 
geometrical  limits  and  the  type  of  joint  between  stationary  and  movable  parts  of  the 
TW.  (This  process  is  similar  to  biological  optimization  of  animal  bones,  which  have 
non-uniform  density  distribution;  it  is  straightforward  that  the  bones  of  birds  have 
the  highest  weight  perfection). 

3.  A  new  dimensionless  parameter  is  introduced  -  load-carrying  factor 
coefficient.  It  is  determined  as  a  ratio  of  LCF  to  the  product  of  characteristic  length 
and  characteristic  load  of  the  structure.  For  things  characteristic  length  corresponds 
to  square  root  of  area  and  characteristic  load  corresponds  to  lift.  This  new  criteria 
allows  for  comparison  of  weight  efficiency  of  different  structures,  and  also  allows  to 
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calculate  absolute  and  fraction  mass  of  the  structure  for  different  geometry  and 
loading.  All  the  necessary  equations  are  obtained. 

Aerodynamic  research  play  a  significant  role  in  the  present  work.  They 
provide  solution  to  two  problems: 

1 .  Identification  of  span  load  distribution  for  TW 

2.  Estimation  of  maximum  lift-to-drag  ratio. 

Because  the  object  of  study  is  unconventional,  three  different  numerical 
methods  were  used  for  aerodynamic  research  in  order  to  increase  assurance  of  the 
results: 

1 .  Linear  bound  vortex  method  (steady  inviscid  incompressible  flow) 

2.  Non-linear  vortex  lattice  method  (nonsteady  inviscid  incompressible  flow) 

3.  Finite  volume  method  (steady  viscious  incompressible  flow). 

Besides  three  different  mockups  were  designed,  built  and  tested  in  the  wind 
tunnel.  In  order  to  provide  validation  and  better  understanding  of  the  results  all 
weight  and  aerodynamic  research  for  TW  were  accompanied  by  the  same  kind  of 
research  for  geometrically  equivalent  trapezoidal  wings. 

Academic  versions  of  MSC  NASTRAN  (license  was  paid  from  financing  of 
the  present  project  REO-1386)  and  STAR-CD  (license  was  already  available  at 
aerohydrodynamics  department)  software  were  used  during  the  research,  together 
with  the  software  developed  in  Samara  State  Aerospace  University  (SSAU). 

The  main  results  of  the  work  are: 

1.  It  was  ascertained  that  telescope  wings  can  have  weight  and  aerodynamic 
efficiency  nearly  equal  to  equivalent  trapezoidal  wings  of  the  same  span  (figures 
9.18,  9.19  and  tables  9.1,  9.2).  With  optimized  geometry  telescope  wings  can  have 
even  better  performance  than  equivalent  trapezoidal  wings. 

2.  The  fundamentals  of  a  new  efficient  approach  to  weight  analysis  were 
developed  for  the  purpose  of  solving  the  stated  problem.  This  approach  is  based  on 
3D-modeling  and  can  be  highly  useful  for  conceptual  design  of  aircraft  of 
unconventional  shapes. 

3.  Project  management  and  financing  through  CRDF  allowed  to  acquire  all  the 
necessary  hardware  and  software  and  to  bring  together  an  efficient  team  of  experts 
capable  to  solve  complex  interdisciplinary  problems. 

The  present  report  is  arranged  in  such  a  way  that  each  part  of  it  can  be  read 
separatly.  That  is  why  some  of  the  results  appear  several  times  in  different  parts. 
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3  Introduction 

Morphing  wing  can  provide  sufficient  improvement  of  weight  and 
aerodynamic  performance  of  aircraft  and  can  fulfill  various  additional  requirements, 
i.e.  they  can  be  adaptive.  For  example,  they  can  be  folding  thus  having  variable 
specific  wing  load. 

Morphing  wing  can  has  different  applications  at  future -technology  flighing 
vehicles.  One  can  find  many  papers  devoted  to  their  study  [21,  22]. 

Usually  morphing  wing  is  an  unconventional  structure.  One  of  the  most 
important  problems  in  the  design  of  aircraft  using  morphing  wing  is  sufficiently 
accurate  estimation  of  morphing  wing  weight  parameters  at  early  stages  of  the 
design. 

Different  methods  of  weight  design  [1,  3,  4,  5,  12,  17,  18,  19,  20]  based  on 
statistical  data  and  simplified  physical  principles  of  structure  operation  give  too 
wide  spreading  of  the  results  even  for  conventional  structures,  as  we  will  show  later. 
Such  approaches  are  not  able  to  account  for  some  principal  features  of  morphing 
wing  -  change  of  the  area  of  lifting  surface,  type  of  the  joint  between  moving  parts 
and  etc.  That  is  why  weight  design  of  morphing  wing  requires  development  of  new 
approaches. 

In  the  present  work  these  approaches  are  developed  in  a  quiet  general  form 
and  their  application  for  telescope  wing,  which  is  a  particular  case  of  the  morphing 
wing,  is  demonstrated. 

In  aviation  the  structural  mass  m*  and  especially  the  structural  mass  fraction 
which  is  mst  to  m0  (aircraft  takeoff  mass)  ratio 


m„ 


m 


st 


m 


(3.1) 


determines  aircraft’s  transport  efficiency  and  even  the  possibility  of  designing  the 
aircraft  with  desirable  performance.  It  can  be  easily  seen  from  existence  equation 
[1,  3]  which  is  used  for  takeoff  mass  estimation  at  early  stages  of  the  design.  It  is 
given  below  with  simplified  set  of  masses  as  an  illustration  only: 

m , 


m0  = 


V 


1  —  m 


St 


m 


sys 


m 


f 


m 


(3.2) 


pp 


where  mpi  -  payload  mass  (for  example  passengers  and  baggage); 
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mst  and  m  -  structural  mass  fraction  and  system  mass  faction,  which  can  be 
estimated  on  statistical  basis; 

mf  and  mpp  -  fuel  mass  fraction  and  propulsion  mass  fraction,  which  depend 

mainly  upon  operational  range  and  speed  and  can  be  calculated  with  quiet  a  good 
accuracy. 

If  the  sum  of  all  relative  masses  in  (3.2)  is  greater  than  one,  then  such  an 
aircraft  could  not  exist.  From  this  comes  the  name  of  the  equation. 

For  some  kinds  of  aircraft  payload  mass  makes  only  10%  of  takeoff  mass  and 
structural  mass  -  25^-30%.  So  on  the  one  hand  one -percent  decrease  of  structural 
mass  causes  several -percent  increase  of  payload  mass,  and  on  the  other  hand  a  small 
error  in  estimation  of  mst  yields  a  great  misestimating  of  the  key  design  parameter 

/Ho¬ 
lt  should  be  mentioned  that  the  problem  of  early  stage  structural  mass 
estimation  is  paid  much  attention  in  aircraft  design  theory.  Basics  of  systematical 
research  of  aircraft  construction  strength  and  weight  efficiency  are  given  in  a  work 
of  F.  R.  Shanley  [19].  Many  of  so  called  “weight  formulas”  for  estimation  of  aircraft 
absolute  and  relative  masses  can  be  found  in  [5,  17,  18,  20]. 

Relative  structural  mass  depends  upon  many  parameters:  geometry, 
accelerations,  load  distribution,  size  and  location  of  cutouts,  material  strength  and 
stiffness  characteristics  and  etc. 

To  take  into  account  dependence  of  internal  structural  loads  from  geometry 
and  external  loads  the  beam  theory  is  often  used  as  a  basis  of  weight  formulas.  Such 
parameters  as  engine  location  (at  the  wing  or  at  the  fuselage),  gloved  wing  and  etc. 
are  accounted  by  multiple  correction  factors  usually  derived  from  statistics  of 
available  aircraft. 

Such  approach  doesn’t  guarantee  adequate  precision  of  weight  calculations  in 
case  of  unusual  scheme,  advanced  constructions  or  dramatically  increased 
dimensions  of  an  aircraft  if  the  “square  and  cube”  law  [18,  19]  is  not  considered  in 
weight  formulas. 

Figure  3.1  presents  an  example  of  wing  mass  fracture  calculation  for  two 
aircraft.  One  is  similar  to  "Iliushin  I1-96MT",  and  another,  transport  airplane,  has 
payload  of  250  000  kg  (project  "Ecolifter"). 
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Figure  3.1  -  Results  of  wing  mass  fracture  estimation 


Different  weight  formulas  give  a  significant  spread  in  values.  To  overcome 
this  shortcoming  V.M.  Sheinin  proposed  his  method  of  multiple  calculations  [18],  in 
which  extreme  estimations  are  discarded  and  intermediate  are  averaged.  This 
method  provides  quiet  high  accuracy  but  its  implementation  is  limited  to 
conventional  structures  and  requires  some  skill  for  usage. 

Figure  3.2  presents  another  example  -  different  load-carrying  schemes  for 
delta  wings.  Known  weight  formulas  do  not  allow  in  principle  to  account  all  features 
of  these  designs. 


Figure  3.2  -  Variants  of  structural  schemes  for  the  delta  wings 
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Thus  a  new  approach  is  needed  in  order  to  increase  the  accuracy  in  aircraft 
structure  design. 

Advantages  in  numerical  simulation  of  structures  and  ubiquitous 
implementation  of  finite  element  method  (FEM),  together  with  determination  of 
some  integral  properties  of  elastic  systems  [10,  13]  allow  create  of  more  general 
approach  for  early-stage  mass  estimation  and  load-carrying  scheme  choice. 
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4  Weight  analisys  of  airframe  structures  based  on  a 
new  metric  "load-carrying  factor" 

4.1  Load-carrying  factor 

Relationship  between  results  of  FEM  simulation  of  the  structure  and  its 
required  mass  can  be  expressed  through  the  specific  criterion  -  "load-carrying 
factor"  -  G,  which  accounts  both  for  magnitude  and  length  of  action  of  internal 
loads  in  the  structure  [10,  12,  13,  15]. 

For  frames 

g=£kk,  (4.i) 

i=  1 

where  i  -  rod  number,  Nt  -  force,  /,■  -  length. 

For  shell  structures 

G  =  Y,RfS,,  (4.2) 

i= 1 

where  i  -  element  number,  Rt  -  equivalent  force  flow,  St  -  element  area. 

For  3D  structures 

G  =  \oeqvdV,  (4.3) 

v 

where  a  eqv  -  equivalent  stress,  V-  structural  material  volume. 

Load-carrying  factor  has  some  interesting  properties  which  can  be  useful  for 
structure  weight  estimation. 

G  is  determined  by  load-carrying  scheme,  i.e.  element  type,  number  of 
elements,  types  of  their  joints  and  arrangement.  For  statically  determinate  structures 
the  value  of  G  is  independent  of  stiffness  ratio  of  load-carrying  elements,  and  for 
statically  undeterminate,  according  to  numerical  simulations,  is  weekly  dependent. 
For  example,  cross-section  areas  of  some  elements  and  their  internal  loads  can 
change  their  order  of  magnitude  during  optimization,  which  result  only  in  few 
percents  of  changing  of  G  value  for  the  whole  structure. 

Theoretically  required  structural  material  volume  VT  for  fully-stressed 
structure  with  definite  load-carrying  scheme  can  be  estimated  through  G  and 
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working  stress  value  after  single  calculation  of  internal  loads  for  some  rational 
initial  distribution  of  element  stiffness  according  to  simple  relationship 


VT 


(4.4) 


(4.5) 


Really,  for  fully-stressed  frame  we  have 

n  In  I  n 

Mb  J  M 

where  F,  -  rod  cross-section  area. 

Similar  result  can  be  obtained  for  shell  structures  using  equivalent  force  flow 
as  a  measure  of  internal  loads  [13]. 

G  value  can  be  calculated  with  quiet  high  accuracy  with  relatively  "coarse" 
FEM  models. 

Using  Vt,  material  density  p  and  full  mass  coefficient  cp,  which  accounts  for 
mass  increase  from  joints  and  non-carrying  elements,  deviations  from  optimal 
material  distribution  in  favor  of  simplicity  and  fabricability  of  the  structure  and  etc., 
one  can  calculate  real  (practical)  structural  mass 

G 


mst  =(p  •  P  =(p  •  p 


[a] 


(4.6) 


or 

q 

mst=<$  — ,  (4.7) 

a 

where  a  -  material  strength-to-weight  ratio,  [o  ]  -  material  allowable  strength. 

In  weight  formula  (4.7)  all  of  the  three  multipliers  determine  nearly 
independent  properties  of  the  structure: 

G  -  geometry,  load-carrying  scheme  and  loads; 
o  -  material  strength; 

(p  -  structural  and  technological  perfection. 

Relationship  (4.7)  can  be  used  for  derivation  of  new  weight  formulas  for  the 
whole  structure  or  its  parts.  In  this  case  cT  will  account  for  material  choice,  G  will 
account  for  all  features  of  load-carrying  scheme  of  the  future  structure,  and  (p 
coefficient  can  be  defined  from  analysis  of  some  already  built  aircraft,  taken  as  a 
prototype 
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(p  = 


mst  G 


(4.8) 


where  sign  marks  values  for  the  prototype. 

If  the  prototype  has  high  weight  perfection  and  no  major  changes  in 
engineering  solutions  (cp)  and  material  choice  (cT )  are  expected,  for  example,  during 
modifications,  then  for  mass  estimation  of  a  new  structure  formula  (4.7)  can  be 
written  with  (4.8)  especially  simple 


mst  = 


m 


St 


(4.9) 


Nowadays  FEM  models  both  for  the  prototype  and  for  a  new  structure  are 
almost  always  available,  so  calculation  of  G  doesn’t  make  a  big  difficulty. 


4.2  Dimensionless  criterion  of  structure 
load-carrying  perfection 


The  disadvantage  of  the  “load-carrying  factor”  criterion  is  its  dimensionality. 
It  makes  difficult  the  usage  of  this  criterion  for  comparison  of  different  engineering 
solutions  at  very  early  design  stages,  when  m0 ,  loads  and  dimensions  of  the  structure 
are  unknown,  thus  it  is  impossible  to  explicitly  estimate  of  structure  weight  fraction. 

While  for  the  broad  range  of  structures,  which  allow  adequate  representation 
by  rod  and  membrane  finite  elements,  it  can  be  shown  that  for  geometrically  similar 
structures  with  the  same  external  nodal  forces  (loads)  the  load-carrying  factor  is 
proportional  to  the  ratio  of  geometrical  similitude  [12].  Thus  the  load-carrying  factor 
is  proportional  to  the  characteristic  dimension  (length)  of  the  structure. 


And  if  we  change  all  nodal  forces  applied  to  the  structure  by  factor  of  k ,  then 
also  the  internal  loads  within  the  structure  will  be  changed  by  factor  of  k,  so  as  its 
load-carrying  factor.  With  provision  for  all  stated  the  above  we  can  introduce  the 
load-carrying  factor  coefficient 


G 

PL  ’ 


(4.10) 


where  G  is  load-carrying  factor  for  the  structure  with  some  specific  load-carrying 
scheme,  dimensions  and  load  distribution;  P  -  characteristic  load;  L  -  characteristic 
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dimension  (length).  Because  the  dimension  of  load-carrying  factor  is  the  product  of 
force  and  length,  the  CK  is  dimensionless. 


4.3  Examples  of  calculations 
of  load-carrying  factor  coefficient  Ck 


Formula  (4.10)  has  a  simple  physical  sense:  CK  is  a  ratio  of  load-carrying 
factor  G  and  characteristic  load  and  characteristic  length. 

Consider  the  basic  problem  (figure  4.1a)  -  transfer  of  force  Pa  by  distance 

ab=l. 

Rod  load  Nah  =  Pa,  load-carrying  factor  for  this  one -rod  structure: 


G  =  Nab-l  =  Pa-l 

and 

(4.11) 

II 

ll 

O 

(4.12) 

Here  force  Pa  is  taken  as  characteristic  load,  and  distance  of  its  transfer  is 
taken  as  reference  length  scale. 

For  the  next  problem  the  direction  of  Pa  is  rotated  by  90°.  A  structure 
consisting  of  two  rods  ac  and  ad  (figure  4.1b)  is  proposed.  With  a  =  45°  CK  =  2,00. 
The  reference  length  scale  is  the  distance  between  point  a  and  rigid  support  -  ab. 


C*-  1,00 


Figure  4.1-  Examples  of  simple  structures 
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Figure  4.1c  show  Mitchell  frame  with  geometrical  limits  ab/cd=2,  for  which 
CK  =  3,41. 

For  a  beam  of  ideal  H-section  with  t/b  =  0,1  (figure  4.  Id)  load-carrying  factor 
computed  from  loads  in  caps  of  spar  only  is  CK  =  10,00. 

Figure  4.2  shows  a  simple  FEM  model  of  three -layer  plate  as  an  example  of 
model  wing  with  planform  be  and  uniformly  distributed  load  p.  The  plate  is  rested 
vertically  at  symmetry  axis.  FEM  model  consists  of  shell  elements  (skin  panels  and 
walls). 


Figure  4.2  -  FEM  model  of  a  three-layer  plate 


For  b  =  c,  which  corresponds  to  aspect  ratio  A  =  1,  and  t=  0,1c  load-carrying 
factor  coefficient  is  CK  =  1,26. 

Figure  4.3  presents  variation  of  CK  with  b/c  =  1;  2;  4;  8;  12,  which  make 
aspect  ratio  A  =  1;  2;  4;  8;  12. 

Figure  4.4  presents  variation  of  CK  with  sweepback  angle  for  aspect  ratio  A  =  8. 
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CK 


Aspect  ratio  A 

Figure  4.3  -  Variation  of  coefficient  Ck  with  wing  aspect  ratio 
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Figure  4.4  -  Variation  of  coefficient  Ck  with  wing  sweepback  angle 


Quite  obvious  relationships  shown  at  figures  4.3,  4.4  illustrate  the  possibility 
of  the  coefficient  CK  to  account  for  geometrical  features  of  the  wing.  Naturally  that 
for  high  aspect  ratio  they  are  in  good  agreement  with  wing  weight  formulas  obtained 
from  beam  theory. 

Dimensionless  load-carrying  factor  coefficient  has  a  definite  physical 
meaning.  Its  value  shows  the  ratio  between  theoretical  mass  of  the  considered 
structure  and  the  mass  of  the  simplest  structure  consisting  of  the  single  rod, 
transferring  the  concentrated  force  along  its  direction. 

The  numerical  value  of  the  load-carrying  factor  coefficient  depends  upon  the 
choice  of  the  structure  characteristic  length  and  load  distribution.  That  is  why  some 
conventions  about  the  choice  of  characteristic  length  and  loading  are  required. 
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4.4  Wing  weight  equations  based  on 
dimensionless  load-carrying  factor  coefficient 


Choosing  \[S  as  a  characteristic  length  of  the  wing  and  accounting  for  the 
above  mentioned  relations  between  load-carrying  factor  and  structure  dimensions 
and  loads,  for  geometrically  similar  wings  we  can  write: 

G  =  CK  -n  ■m0  ■  g  -J~S,  (4.13) 

where  n  -  load  factor,  S  -  wing  area. 

Relation  (4.13)  can  be  used  for  determination  of  CK  from  the  calculation  for 
some  geometrically  similar  wing  (values  for  this  wing  are  marked  with  "  "  sign). 


G* 

*  *  /  0* 

n  ■m0  ■  g-^S 


(4.14) 


In  (4.13)  and  (4.14)  4s  acts  as  a  dimensional  coefficient  of  geometric 
similarity  -  characteristic  length,  and  product  n-mQ  ■  g  determines  the  load. 
Dimensionless  magnitude  CK  accounts  for  wing  shape,  structure  and  load 
distribution. 

Load-carrying  factor  coefficient  CK  is  equal  to  load-carrying  factor  for  the 
wing  of  some  definite  shape  and  structure,  with  unit  area  and  unit  net  loading. 
Accounting  for 

m 

S  =  -2-,  (4.15) 

Po 

where  pa  -  wing  specific  load, 

We  can  write  (4.13)  as 


G  =  CK  n  ■ m0-g ■ 


(4.16) 


Then  from  (4.7)  wing  mass  fracture  can  be  found  as 


m 


wing 


9  r 

=  ~CK-n  g- 

a 


(4.17) 


or 
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mwing  =~=CK  n'  (4-18) 

Absolute  mass  of  the  wing  structure  will  be 

mmng=^:CK  -n-m0  -g-Vs.  (4.19) 

Thus,  load-carrying  factor  coefficient  is  a  dimensionless  value,  which  allows 
calculation  of  structure  absolute  mass  and  mass  fracture  at  early  stages  of  design 
and,  what  is  the  most  important,  it  allows  comparison  of  different  load-carrying 
schemes. 


4.5  Load-carrying  factor  of  reinforced  panels 


On  modem  aircraft  stmctures  (especially  in  high-loaded  wing  structures)  the 
main  load-carrying  element  is  stringered  skin.  The  most  important  parameters  for 
such  stmctures  are  orientation  of  the  stringers  and  relative  skin  thickness  8: 


8  = 


^  skin 


reduct 


(4.20) 


where  8Skin  -  skin  thickness;  &reduct  =  8  skin  +Fs  fbs  -  reduced  panel  thickness;  Fs  and 
bs  -  cross-section  area  and  pitch  of  the  stringers. 

Optimal  value  of  8  rising  from  stability  requirements  lies  in  a  narrow  range 
from  0.4  to  0.6,  where  0.4  is  for  panels  with  high  compressive  loads  acting  parallel 
to  the  stringers. 

According  to  load-carrying  factor  definition  (4.1,  4.2,  4.3)  and  its  further  use 
for  weight  estimation  (4.7),  two  strategies  for  load-carrying  factor  calculation  can  be 
proposed  for  stringered  panels. 

Strategy  I.  Stringered  panel  is  considered  as  substructure,  consisting  of  the 
shell  skin  and  a  number  of  rods.  Then: 


Gpan=Gstln  +  Gr  =o  Z -V,Ua  +KK,  (4.21) 

where  Gskin  -  skin  load-carrying  factor,  G>  -  rods  load-carrying  factor,  Vskin  -  skin 
material  volume,  Vr  -  rods  material  volume,  a  eJ?n  -  skin  equivalent  stress,  of.  -  rods 
equivalent  stress. 

In  fact,  this  strategy  for  a  panel  realizes  the  most  general  expression  (4.3). 

Strategy  II.  The  alternative  is  the  calculation  of  panel  load-carrying  factor 
using  the  values  of  maximum  equivalent  stress  in  the  panel  and  of  its  full  volume: 
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G  =<5eqv 

pan  max 


v 


pan  ’ 


(4.22) 


where  o  ^'x  -  maximum  equivalent  stress  in  the  panel,  Vpan=  Vskin  +  Vr  -  panel 
volume. 

The  purpose  of  calculating  panel  load-carrying  factor  in  this  way  is  the 
estimation  of  reinforced  panel  required  mass  from  strength  requirements  with 
prescribed  value  of  8  : 


R 


eqv 


1 


X2  +  — 


xy  „r 

^  +  3  — 


9  reduct 


(4.23) 


[a  ]  [o  ] 

where  R  eqv  -  equivalent  force  flow  in  a  reinforced  panel,  X,  Y,  T -  normal  and  shear 
force  fluxes  in  the  panel. 

Notice  that  for  such  choice  of  panel  reduced  thickness  5reduct  rod  elements  can 
be  underloaded. 

Special  research  has  shown  that  wing  load-carrying  factor  is  highly  dependent 
on  orientation  of  stringers.  Calculation  of  load-carrying  factor  G  according  to 
strategy  I  is  rational  in  case  when  relative  skin  thickness  8  or  stringer  orientation 
can  be  change  during  the  process  of  mass  optimization.  Usage  of  strategy  II  is 
rational  for  weight  estimation  of  structures  with  different  load-carrying  structures, 
defined  by  stringer  orientation.  Implementation  of  both  strategies  allows  to  estimate 
possibilities  for  load-carrying  structure  optimization. 


4.6  Application  of  load-carrying  factor  coefficient 
for  estimation  of  weight  efficiency 
of  wing  load-carrying  structure 

Consider  the  application  of  load-carrying  factor  coefficient  for  comparative 
analysis  of  load-carrying  structures  of  delta  wings  as  example.  For  this  purpose 
structures  d,  e,  b  from  figure  3.2  are  chosen,  they  are  numbered  1,  2,  3  respectively. 
Structures  1  and  2  are  widely  used  and  are  known  as  concurrent. 

In  the  first  structure  spars  are  span-oriented,  which  gives  the  shortest  paths  of 
load  transfer  in  the  structure,  but  the  longest  spars  in  this  case  are  located  in  the  area 
of  low  structural  height. 
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The  second  structure  with  spars  and  stringers  coming  from  the  wing  tip  has  its 
own  load-carrying  and  technological  advantages.  This  structure  provides  load 
transfer  through  the  areas  with  high  structural  height,  cap  of  spar  bevel  is  constant 
along  the  span.  The  disadvantage  is  structurally  unavoidable  otherstrengh  of  wing 
panels  near  the  tip.  Root  rib  for  this  structure  should  be  more  thick  than  for  the  first 
one  because  of  the  kink  of  longitudinal  load-carrying  elements. 

Without  quantitative  analysis  it  is  impossible  to  choose  the  best  of  these  two 
structures,  neither  it  is  possible  to  make  an  accurate  estimate  for  absolute  mass  and 
mass  fracture  for  these  wings. 

The  third  structure  uses  only  one  spar  for  torque  Mx  transfer  (from  layout 
requirement).  It  is  interesting,  how  inefficient  this  structure  in  comparison  with  the 
first  two. 

For  calculations  of  load-carrying  factor  coefficient  is  such  problems  it  is 
rational  to  use  advanced  implementation  of  FEM  (MSC  NASTRAN,  ANSYS, 
ARGON,  RIPAK  and  etc.). 

Figure  4.5  presents  FEM  models  for  these  wings.  Airfoil  is  the  same  for  all 
three  wings,  lenticular  with  relative  thickness  T  =  10%.  Skin  panel  and  walls  are 
modeled  with  isotropic  elements.  Depending  upon  analysis  strategy  and  features  of 
load-carrying  structures  rod  elements  are  added,  to  model  caps  of  spar  and  stringers. 

Aerodynamic  loading  is  considered  to  be  uniform,  because  the  main  purpose 
of  the  wing  load-carrying  scheme  is  transfer  of  aerodynamic  loading  to  the  fuselage 
and  for  correct  comparison  of  load-carrying  structures  some  standardization  of  load 
distribution  is  needed. 

Geometry  is  the  same  for  all  wings:  wing  span  b  is  equal  to  root  chord  c  , 
belly  section  -  0,222  b. 

Results  are  presented  in  table  4. 1 . 
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Wing  1 


Wing  3 


Figure  4.5  -  FEM  models  of  the  wings 


Table  4.1  -  Load-carrying  factor  coefficients  for  the  wings 


Wing 

Membrane 

structures 

Panel  structures 

Strategy  I 

Strategy  II 

5  =  0,6 

5  =  0,5 

5  =  0,4 

5  =  0,6 

5  =  0,5 

5  =  0,4 

1 

1,62 

1,68 

1,70 

1,71 

1,84 

1,94 

2,07 

2 

1,68 

1,76 

1,78 

1,81 

1,83 

1,89 

1,98 

3 

2,55 

2,69 

2,75 

2,83 

2,68 

3,03 

3,56 

Models  composed  of  membrane  elements  only  were  used  for  estimation  of 
minimal  value  of  load-carrying  factor  coefficient  which  can  be  achieved  for 
prescribed  constraints  conditions.  As  expected,  values  for  wings  1  and  2  are  very 
close.  Weight  efficiency  of  the  wing  3  is  57  %  less  (rod  elements  were  used  to 
model  the  central  beam). 

Calculations  were  performed  in  a  following  way.  From  FEM  model  for  each 
membrane  element  of  the  wing  skin  its  area  Sm  and  volume  Vm  =  Sm  ■  bskin  were 
defined.  Then  volume  of  rods  for  the  element  was  calculated 

vr=vm~.  (4.24) 
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Figure  4.6  -  Modeling  of  stringered  panel  by  shell  and  rods 


From  known  lengths  a  and  c  (figure  4.6)  cross-section  areas  for  the  rods  were 

found: 


Fr=-^—.  (4.25) 

/  +/  v  ' 

a  c 

On  nodes  1-2  and  3-4  (in  the  direction  close  to  stringers)  rod  elements  with 
cross-section  area  Fr  and  lengths  la  and  lc  were  introduced.  Then  static  calculations 
were  performed,  following  by  calculations  of  load-carrying  factor  G  according  to 
strategies  I  and  II.  Load-carrying  factor  coefficient  was  calculated  from  (4.10), 
where  sum  of  all  nodal  forces  (loads)  in  FRM  model  were  taken  as  a  characteristic 
load: 


Ck  — 


£/> -  Vs 


(4.26) 


j 

where  Pj  -  nodal  forces,  applied  to  the  console  of  wing;  S  -  area  of  wing  together 
with  belly  section. 

From  table  4.1  it  follows  that  the  results  obtained  from  strategy  I  are  quiet 
close  for  all  8  and  there  is  only  slight  difference  from  membrane  models,  which 
give  the  estimation  for  minimal  possible  value  of  CK  for  prescribed  conditions  of 
wing- fuselage  joint.  Results  obtained  from  strategy  II  are  adequate  in  accounting  for 
specific  of  stringered  panels  straining,  thus  providing  the  possibility  for  objective 
comparison  of  load-carrying  structures  with  different  values  of  8  and  calculation  of 
absolute  mass  and  mass  fracture  for  such  wings  from  formulas  (4.19,  4.18). 
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5  Using  3D  FEM  modeling 
for  weight  calculations 

5.1  Theoretical  background 

Weight  analysis  of  morphing  wing  requires  consideration  of  unusual  external 
shapes,  for  which  it’s  very  difficult  to  propose  rational  airframe  structures  and  create 
FEM  models  like  those  that  shown  at  figure  3.2. 

That  is  why  for  weight  analysis  of  morphing  wing  it  is  reasonable  to  use  3D 
continuum  FEM  models  which  fill  all  allowed  volume  with  the  account  for 
geometrical  constraints,  i.e.  they  have  the  external  shape  as  that  of  the  wing  and  do 
not  fill  hollows  in  which  it  is  not  allowed  to  place  load-carrying  elements  due  to 
some  reasons.  Such  models  contain  all  possible  airframe  structures,  including  the 
optimal  ones. 

Probably  for  the  first  time  the  idea  of  this  aproach  was  stated  in  [10].  Basic 
algorithm  for  density  distribution  optimisation  for  3D-models  presented  in  [14]. 

The  theoretically  required  mass  of  the  optimal  structure,  contained  within  3D- 
model,  can  be  estimated  as  following. 

Consider  some  hypothetical  material  with  non-uniform  density  distribution. 
Consider  its  strength  -  allow  stress  [o],  and  coefficient  of  elasticity  E  are  linear 
functions  of  density: 

[a]=p[a],  (5.1) 

E=pE,  (5.2) 

where  [<T ]  and  E  are  strength  and  coefficient  of  elasticity  for  unit  density. 

Assign  some  amount  of  this  material  m  for  creation  of  3D-model.  Mesh  it  with 
finite  elements  of  “solid”  type  and  assign  the  same  density  and  coefficient  of 
elasticity  for  all  of  the  elements: 
m 

Po,=-,  (5-3) 

£<»=PoA  (5-4) 

where  subscripts  denote:  0  -  initial  distribution  ,  i  -  element  number,  V  -  total 
volume  of  all  finite  elements  of  the  3D -model. 
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Auxiliary  optimization  problem  can  be  formulated  as  following: 

Minimize  potential  strain  energy  of  3D-model 

U  — >  min  (5.5) 

under  the  constraint  of  its  constant  mass 

m  =  const.  (5.6) 

with  the  assumption  that  densities  of  elements  can  vary,  i.e.  they  are  considered  as 
design  variables: 

X;-  P;  ?  i  1  ,...n  ,  (5.7) 

where  n  is  the  number  of  elements  if  the  3D -model. 

To  solve  this  problem  we  use  the  main  idea  of  [10]. 

Constrain  the  variables  p  ■  with  mass  conservation  (5.6) 


£  px=«. 


(5.8) 


i= 1 


where  V,  -  volume  of  element. 

Compute  the  initial  uniform  distribution  of  the  filler  density 

(  n  V* 1 


Poi  =  m  ' 


Lr, 


1=1 


(5.9) 


where  V,  is  the  volume  of  the  element. 

Define  the  principal  stresses  within  elements  G10i,  G  2()i ,  G  30/  and  write  down 
the  relationship  for  3D -model  energy 


y  =  E 


n  (j 


eqv  o  i 


V, 


(5.10) 


N2  PiE 

Relation  (5.10)  assumes  that  stresses  in  elements  are  considered  as  constant, 
and  densities  p;  can  vary. 

To  find  the  optimal  values  of  densities  under  proposed  assumptions  we  use 
Langrangian  multiplier  method  [11]  and  write  the  following  system  of  equations 


dU  .  dcp  ,.  , 

—  +  ^  —  =  0,(z  =  1,2,...,  n) 

dp*  3P/ 

n 

<P(P/)  =  Epf^ -m  =  0, 

i= 1 


(5.11) 


where  X  -  Langrangian  multiplier  ,  (p  -  auxiliary  function,  expressing  the  mass 
conservation  (5.8). 

Differentiate  with  the  account  for  (5.10)  and  obtain 
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eqvoi  i 


2Epf 

£p^-™  =  o 

L/'=l 

Whence  it  follows  that 

a. 


fXF  =0,(1  =  1,2 


P,  = 


eqvoi 

42EX 


(5.12) 


(5.13) 


and  with  the  account  for  (5.8) 

P  if  =-^— ^ - %„•.  (5-14) 

Vo  U 

/  v  egvoz  z 
1=1 

where  pl(.  is  the  found  new  density  distribution,  instead  of  some  initial  po(  set,  for 
example,  by  (5.3). 

With  stresses  a  eqvfn  corresponding  to  elements  density  distribution  po/  and 
new  density  distribution  pl(. ,  obtained  from  solving  (5.11)  the  extremum  problem 
(5. 5 -5. 7)  potential  strain  energy  of  3D -model  will  be  lower  and  become  Ul_0<U0, 
where  U0  -  potential  energy  of  the  initial  3D -model  with  stresses  g(123)0j.  and 
elements  density  poi ;  Ux_0  -  energy  of  3D -model  with  stresses  a  ;  and  new 
density  distribution  p  u .  Energy  UX_Q  can  be  calculated  by  substitution  of  p  u  into 


(5.10). 

But  3D -model  with  a  new  density  distribution  pl(  will  have  another  stress 
distribution  G(l  2  3)l/ ,  which  will  satisfy  equilibrium  and  deformation  consistency 
conditions.  In  a  new  3D -model  stress  distribution  c  (l  2  3)  oi  will  satisfy  equilibrium 

condition  only,  but  not  the  deformation  consistency.  Actual  potential  strain  energy 
of  3D -model  with  new  distribution  of  density  pb  will  be  Ul  and  according  to  the 
principle  of  minimal  strain  energy 

UX<U^.  (5.15) 

Using  the  new  distribution  pl;.  as  the  initial  one  and  making  successive 

calculations  of  structure  stress  state  according  to  (5.14),  we  get  a  sequence  of 
structures  with  energies 

U0>Ul>U2>...  (5.16) 
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This  monotonically  decreasing  sequence  is  limited  from  below,  because  the 
strain  energy  is  always  positive,  and  thus  it  is  convergent.  Reduction  of  the  energy  to 
the  minimal  possible  value  means  a  search  for  density  distribution  into  3D -model 
that  gives  the  maximal  structure  stiffness. 

We  can  also  prove  [11]  that  this  iteration  process  yields  to  the  structure  with 

n 

minimal  value  Voet?v/-  Vi ,  i.e.  with  minimal  load-carrying  factor  G,  calculated  by 

i= 1 

(4.3). 

The  knowledge  of  this  value  allows  as  to  determine  from  (4.7)  the  minimal 
theoretically  required  mass  of  the  structure,  which  can  be  designed  with  the  account 
for  geometrical  constraints. 

In  the  proposed  iteration  process  high-loaded  elements  will  increase  their 
density  and  low-loaded  will  decrease.  And  variation  of  density  will  continue  until  all 
the  elements  with  non-zero  density  will  have  will  have  equal  specific  strain  energy 
to  density  ratio.  Density  of  some  elements  will  be  nearly  zero,  which  indicate  non¬ 
rationality  of  placement  of  load-carrying  elements  in  these  zones. 


5.2  Methodology  of  testing 

During  the  calculations  (5.14)  new  values  of  element  densities  are  assigned 
proportionally  to  the  equivalent  stresses  within  them  Because  for  weight  analysis 
we  are  interested  in  the  minimal  possible  value  of  the  load-carrying  factor, 
corresponding  to  the  given  density  distribution,  these  densities  can  be  evaluated  with 
more  simple  relation: 

^  6([V  oi  /Cl  ^7\ 

P*  =  TT  (  } 

As  an  example  consider  the  calculation  of  load-carrying  factor  for  3D -model 
of  the  wing  structure  (figure  5.1). 
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Figure  5.1  shows  wing  box  structure  with  b/c  =  8,  t/c  =  0. 1,  fixed  in  the  root 
section  and  loaded  with  uniformly  distributed  load.  Calculation  of  theoretically 
optimal  thin-wall  structure,  consisting  of  isotropic  skins  and  shear-carrying  walls, 
yields  for  the  given  geometrical  constraints  the  value  of  load-carrying  factor 
coefficient  CK  =  19,58  (figure  5.2).  Calculation  for  the  wing  model  consisting  of 
solid  elements  in  case  of  uniform  density  distribution  yields  CK  =  28,10.  But  in  the 
iteration  optimization  process  for  density  distribution  after  only  10  iterations  we 
obtain  CK  =  22,1 1  with  8-layer  meshing  of  3D  wing  model,  which  can  be  considered 
reasonable  for  estimation  purposes.  As  it  was  expected,  in  course  of  optimization 
density  of  the  high-loaded  surface  layers  was  increased. 
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c  K 


wing  box  o  solid  model  Aspect  ratio  A 

Figure  5.2  -  Load-carrying  factor  coefficient  of  wing  for  different  aspect  ratio 
Investigation  for  aspect  ratio  A  =  2,  4,  8,  12  is  shown  at  figures  5.2  and  5.3. 
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Figure  5.3  -  Comparison  of  load -carrying  factor  coefficient  calculations  for  thin -wall 
structure  and  for  3D-solid  model  with  variable  density 
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From  8 -layer  model  the  value  of  CK  is  12%  higher  comparing  with  ideal 
structure,  as  it  was  expected.  It  is  explained  by  the  fact  that  for  the  optimal  density 
distribution  in  8 -layer  3D  model  the  real  wing  thickness  will  be  7/8  t. 

As  the  number  of  layers  in  the  model  increased,  the  results  of  CK  calculation 
with  3D  model  approach  asymptotically  to  the  theoretical  optimum.  But  because  of 
the  necessity  to  consider  multiple  variants  of  the  structure,  only  8-layer  models  were 
used  in  all  calculations  of  this  research. 

Thus  3D-modelling  allows  to  obtain  the  estimation  of  theoretically  necessary 
mass  of  the  structure  with  reasonable  accuracy  and  considerably  simplifies 
preparation  of  FEM  models. 
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6  Full  weight  coefficient  estimation 


Full  mass  coefficient  was  calculated  for  the  wing  of  maneuvering  aircraft 
Yakovlev  Yak- 130.  General  drawing  of  the  aircraft  is  shown  at  figure  6.1.  Finite- 
element  model  of  the  wing  developed  for  the  calculation  of  (p  is  shown  at  figure  6.2. 
Calculations  provided  the  value  (p  =3,3. 


Figure  6.1-  Layout  of  maneuver  airplane  Yak- 130 
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Figure  6.2  -  Finite-element  model  of  Yak-130  wing 

Because  there  are  no  prototypes  of  telescope  wings  available  nowadays,  now 
the  full  mass  coefficient  is  being  calculated  for  the  flap  of  a  liner  wing.  It  was 
chosen  as  one  of  the  possible  prototypes  because  its  design  includes  rails  and  contact 
pairs.  3D-model  for  calculation  of  cp  is  shown  at  figure  6.3. 

For  this  structure  full  mass  coefficient  about  5  was  obtained. 

Calculation  for  flap  was  performed  with  one  load  case  only.  This  of  course 
yields  underestimated  value  of  load-carrying  factor  for  this  structure,  especially  for 
the  rail,  and  consequently  results  in  overestimated  value  of  full  mass  coefficient. 

The  same  result  was  obtained  for  slat  (figure  6.4). 

The  method  for  calculation  of  full  mass  coefficient  for  movable  aircraft 
structures  require  additional  investigation. 


REO-1386  Research  of  Morphing  Wing  Efficiency 

Final  Technical  Report  /  Samara  State  Aerospace  University,  AVIKON  Institute 
1  October  2003  -  1  June  2004 


35 


Figure  6.3  -  Finite-element  model  of  a  flap 


Figure  6.4  -  Finite-element  model  of  a  slat 
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7  Computational  research  telescope  wings 

aerodynamics 


7.1  Brief  description  of  computational  method  of 
aerodynamic  calculations 


The  method  for  calculations  of  telescope  wing  aerodynamics  basing  on 
discrete  vortex  method  with  linear  formulation  [2,  6,  8]  was  used  for  calculations  of 
angle  of  attack  derivatives  of  lift  coefficient  and  pitching  moment  coefficient,  and 
drag  due  to  lift  factor.  Zero  lift  drag  was  calculated  with  empirical  formula,  basing 
upon  equivalent  flat  plate  skin  friction. 

Main  ideas  of  discrete  vortex  method  are  widely  known  and  used  for 
aerodynamic  calculations.  The  wing  is  replaced  by  N  chordwise  strips,  and  each 
strip  is  divided  by  n  panels.  At  each  panel  an  oblique  horseshoe  vortex  and  an 
impermeability  point  (collocation  point)  are  placed.  Computational  scheme  "1/4  — 
3/4"  was  used  in  the  present  method.  In  this  scheme  the  attached  horseshoe  vortex  is 
located  1/4  of  the  panel  chord  from  the  leading  edge  of  the  panel,  and  collocation 
point  3/4  of  the  panel  chord  from  the  leading  edge  of  the  panel  (figure  7.1). 


z 


X 


Figure  7. 1  -  Computational  scheme  of  discrete  vortex  method 


Impermeability  condition  satisfied  in  the  collocation  points  yields  to  the 
system  of  linear  algebraic  equations  about  unknown  intensities  of  the  vortices: 


REO-1386  Research  of  Morphing  Wing  Efficiency 

Final  Technical  Report  /  Samara  State  Aerospace  University,  AVIKON  Institute 
1  October  2003  -  1  June  2004 


37 


m  _ 

Y,(wylJ+AwylJ)T?=-2K,  j  =  \,m ;  (7.1) 

1=1 

where  wyiJ- =wyi^0ij,q0ijA[l)  is  dimensionless  velocity  induced  by  vortex  with 

number *  1  in  the  collocation  point  with  number  /; 

Awyy  =  wyi^0ij,  Aq0iJ>  A^)  is  additional  dimensionless  velocity  induced  in  the 

collocation  point  j  by  the  mirrored  vortex  symmetrical  to  vortex  i; 

^0ij,q0ij  are  coordinates  of  collocation  point  in  the  coordinate  system  with  the 

origin  in  the  middle  of  the  vortex; 

is  sweepback  angle  of  the  vortex. 

Solving  system  (7.1)  we  obtain  T“  and  the  integral  aerodynamic  parameters 
of  the  wing  can  be  found  as  following: 


A  rn 

(7.2) 

A  171  r 

f  _  >  ra 

^  i  r  ’ 

A  i= i  CA 

(7.3) 

-a  C 

“  CL  ’ 

(7.4) 

i  /  n 

c  =  fr 

4sec  N  c  *-1  1  ’ 

iv  Sec  »=1 

(7.5) 

-  Ca  c 

■p  _  ysec  sec 

Cac  ’ 

(7.6) 

where  x.  is  the  coordinate  of  the  middle  point  of  i  -th  horseshoe  bound  vortex. 

Zero  lift  drag  coefficient  of  the  wing  was  calculated  as 

C0q  =  2C7 (0,93  +  2,8  7) (l  +  5  TM4), 

(7.7) 

flat -plate  skin  friction  coefficient  was  calculated  as  the  sum 

2Cf  =  2Cn  —  +  2Cf2—, 

f  fl  s  ft  s 

(7.8) 

flat -plate  skin  friction  coefficient  for  each  part  of  the  wing  was  calculated  according 
to  Reynolds  number  under  the  assumption  of  turbulent  boundary  layer 
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2C,= 


0,91 


(1  +  0,1M2) 


log 


Re 


\ 


-|  2,58  * 


(l  +  0,lM2) 


yj 


(7.9) 


Reynolds  number  for  each  part  of  the  wing  was  based  on  its  mean  chord 
V  c 

Re;  =  - — . 
v 

Coefficient  of  kinematic  viscosity  was  taken  from  standard  atmosphere  for  the 
altitude  of  0  m. 

Maximum  lift-to-drag  ratio  was  calculated  as 
1 


^Aviax 


K 


7.2  Testing  of  computational  program  for  calculation  of 
spanwise  circulation  distribution 

In  order  to  validate  calculation  results  obtained  with  discrete  vortex  method 
[2,  6]  a  number  of  test  calculations  were  performed  for  trapezoidal  wings, 
geometrically  similar  to  the  considered  telescope  wing.  Table  .1  presents 
comparison  for  two  wings  of  the  same  aspect  ratio  10  and  same  leading  edge 
sweepback  angle,  obtained  from  linear  discrete  vortex  method  [8]. 


Table  7.1  -  Computed  angle  of  attack  derivative  of  lift  coefficient  for  two  trapezoidal  wings 


Wing 

cL , 

from  paper  [3], 

mesh  1274 

CL  , 

present 

calculations, 

mesh  200?  20 

Relative  error,  % 

A  =  10;  A,  =  1.0;  A/e=  0.17 

4,951608 

4,842948 

2,24 

T  =  10;  k  =  0.5;  A/e=  0.17 

5,069728 

4,983007 

1,74 

Figures  7.2  and  7.3  show  spanwise  distribution  of  dimensionless  circulation. 
A  good  agreement  between  the  results  can  be  seen.  The  difference  in  the  values  of 
CLa  ,  which  is  visible  in  table  7.1,  can  be  explained  by  the  difference  in 

computational  grids.  It  is  obvious  that  the  more  fine  grid  used  in  present  calculations 
provides  more  accurate  result. 
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t  0,2  0,4  0,<5  0,8  z 

• -calculafconal  data,Kalinin,  TSAGI  [3j; 

-  -calculationa!  data 


Figure  7.2  -  Spanwise  load  distribution  for  the  wing  with  unit  taper  ratio 


— calculational  data 


Figure  7.3  -  Spanwise  load  distribution  for  the  wing  with  taper  ratio  0.5 

Figures  7.4  and  7.5  shows  comparison  of  computational  results  with 
experimental  data. 
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*  -  experimental  data,  KaMtx  TSAG!  [3] ; 

--calculational  data 

Figure  7.4  -  Spanwise  load  distribution  for  the  wing  with  taper  ratio  0.4 


*  -  experimental  data,  Kaliniix  TSAGI  [3] ; 
—  -calculatLona]  data 


Figure  7.5  -  Spanwise  load  distribution  for  the  sweptback  wing  with  taper  ratio  0.4 

Comparison  with  other  calculations  and  experimental  data  has  shown  a  good 
agreement. 

Table *  1 2  presents  aerodynamic  performance  for  the  trapezoidal  wing  with 
aspect  ratio  A  =  10,964  and  taper  ratio  X  =  0,5  (wing  W_2b)  and  for  the  telescope 
wing  (wing  W_2a)  calculated  with  discrete  vortex  method.  Zero  lift  drag  coefficient 
was  calculated  for  the  airspeed  35  m/s. 
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Table  7.2  -  Calculated  aerodynamic  performance  of  trapezoidal  and  telescope  wing 


Wing 

Aspect 

ratio 

Area, 

m2 

CL 

—a 

XF 

CD 

vo 

K 

L/D 

/  max 

Mesh 

W_2a 

10.800 

0,012 

5,0127 

-0,0240 

0,2746 

0,0185 

0,0301 

21,179 

96?20+48?20 

W_2b 

10,964 

0,012 

4,9938 

-0,0259 

0,2976 

0,0229 

0,0297 

19,185 

100720 

Calculations  were  performed  for  the  models  which  were  tested  in  the  wind- 
tunnel  (see  section  8).  Table  7.3  and  figure  7.6  present  results  for  four  telescope 
wings  with  different  span  of  the  movable  part. 


Table  7.3  -  Calculated  aerodynamic  performance  of  four  different  telescope  wings 


Wing 

Aspect 

ratio 

Area, 

m2 

CL 

—a 

XF 

CD 

Do 

K 

LID 

/  max 

Mesh 

W_05 

3,273 

0,110 

3,3363 

-0,0141 

0,2437 

0,0124 

0,0979 

14,338 

125720+25720 

W_15 

4,923 

0,130 

3,9739 

-0,0191 

0,2755 

0,0127 

0,0653 

17,357 

125720+75720 

W_20 

5,786 

0,140 

4,1887 

-0.0209 

0,2856 

0,0128 

0,0557 

18,717 

100720+80720 

W_30 

7,562 

0,160 

4,5021 

-0.0236 

0,3002 

0,0130 

0,0427 

21,211 

75720+90720 

Figure  7.6  -  Spanwise  load  distribution  for  four  different  telescope  wings 


The  calculations  have  shown  that  the  increase  of  the  relative  span  of  the 
movable  part  of  the  telescope  wing  results  in  the  more  sharp  “kink”  of  the  load 
distribution  at  the  joint  between  the  two  parts. 
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7.3  Results  of  the  calculations 


Results 


of  the  calculations  for  six  different  wings  shown  at  figure  7.7  are 


considered. 


b/2=12 


Ti¬ 

ll 

Q 


n 

N 


Ti¬ 

ll 

u 


II 


Wing  1  a 


Wing  lb 


Wing  2  a 


Wing  2  b 


Wing  3  a 


Wing  3  b 


Ct  =  0,75 


Ct  -  0. 5 


C,  =  0.25 


Figure  7.7  -  Geometry  of  the  wings 

The  table  7.4  below  shows  aerodynamic  coefficients  for  all  of  the  six  wings. 
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Table  7.4  -  Aerodynamic  coefficients  for  all  of  the  six  wings 


b/2  -  1 8 

Wing 

0,75 

T  =  9.818;C4  =  4.9 197;  if  =  0.0354; 

CDo  =  0.00636;  L/DMAX  =  33.290 

la 

— 

T  =  10.058;Q  =  4.9362;  if  =  0.0323; 

CDo  =  0.00643;  L/Dmax  =34.692 

lb 

0,5 

T  =  10.800;C,  =  5.0183;if  =  0.0301; 

CDo  =0.00639 ;L/Dmax  =36.051 

2a 

— 

T  =  10.941;C,  =  5.0238;  if  =  0.0297; 

CDg  =  0.00653;  L/Dmax  =35.895 

2b 

0,25 

A  =  12.000;  C4  =4.9096;K  =  0.0311; 

CDg  =  0. 00638;  L/  Dmax  =  35.461 

3a 

— 

A  =  12.000;  C4  =5.1 146;^  =  0.0272; 

CDg  =  0.00663;  L/Dmax  =37.257 

3b 

Figure  7.8  shows  the  spanwise  load  distribution  for  the  telescope  wings. 
Figures  7.9  and  7.10  present  the  spanwise  distribution  of  dimensionless  velocity 
circulation  for  wings  1  and  5. 


Figure  7.8  -  Spanwise  load  distribution  for  the  telescope  wing  with  semispan  18 
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Figure  7.9  -  Spanwise  load  distribution  for  trapezoidal  wing  "  1" 


Figure  7. 10  -  Spanwise  load  distribution  for  trapezoidal  wing  "5" 

As  one  can  see  from  the  table  7.4,  the  telescope  wing  with  ct  =  0.5  (wing  2a) 
provides  more  lift  and  has  better  aerodynamic  efficiency. 

Figure  7.1 1  presents  spanwise  load  distribution  for  the  rectangular  wing  with 
aspect  ratio  A  =  6,0. 
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Figure  7.11-  Spanwise  load  distribution  for  the  rectangular  wing  with  aspect  ratio  A  =  6,0 


Table  7.5  presents  calculated  aerodynamic  performance  for  the  rectangular 
wing  with  aspect  ratio  A  =  6,0. 

Table  7.5  -  Calculated  aerodynamic  performance  for  the  rectangular  wing  with  aspect  ratio  A  =  6.0 


Wing 

Aspect 

ratio 

Area,  m2 

MX 

—a 

Xf 

CD 

Do 

K 

Mesh 

W_6 

6,0 

96,0 

4,2179 

-0,0175 

0,2388 

0,0078 

0,0601 

23,112 

192720 

7.4  Results  of  numerical  simulations 
for  the  test  model  with  pressure  orifices 

For  computational  identification  of  loads  and  pressure  distribution  a  CFD 
model  was  created  using  Star-CD™  software.  This  software  is  developed  by  CD- 
Adapco®,  permanent  academic  license  was  provided  to  SSAU  by  Russian 
representation  of  distributing  company  CAD-FEM  GmbH.  This  software  provides 
solution  to  general  governing  equations  of  fluid  dynamics  (Navier-  Stokes  equations) 
with  finite  volume  method. 

CAD  model  of  the  telescope  wing  is  shown  at  figure  7.12.  It  has  the  same 
dimensions  as  the  wind-tunnel  test  model,  designed  for  study  of  pressure 
distribution,  which  is  described  in  the  next  section.  Figure  7.13  shows  the  details  of 
the  surface  mesh  in  the  joint  area.  For  both  parts  of  the  wing  mesh  is  collocated  to 
the  joint  in  order  to  provide  accurate  resolution  of  the  sharp  gradients  of  the  flow 
parameters  which  are  expected  in  this  region.  The  size  of  the  mesh  cell  near  the  joint 
is  about  0,2  mm  in  spanwise  direction. 
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Figure  7. 12  -  Geometrical  model  of  the  telescope  wing 


Figure  7. 13  -  Mesh  in  the  joint  area 
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Figure  7.14  shows  the  general  view  of  the  computational  domain  around  the 
telescope  wing,  and  figures  7.15  and  7.16  the  same  domain  as  it  is  seen  when 
viewed  along  Z,  X  and  Y  axis  of  global  Cartesian  coordinate  system.  Global 
Cartesian  coordinate  system  is  connected  to  the  wing  and  the  directions  of  its  axis 
correspond  to  usually  assumed  in  aerodynamics  -  X  is  pointing  backwards 
(streamwise)  along  the  root  chord,  Y  is  pointing  up,  Z  is  pointing  in  spanwise 
direction. 

Only  one  half  of  the  wing  was  considered  in  order  to  reduce  computational 
expenses,  i.e.  the  flow  was  considered  to  be  symmetrical.  Distances  between  the 
wing  and  the  boundaries  of  the  computational  domain  were  taken  according  to 
standard  practice  in  CFD  simulations.  Inlet,  upper  and  lower  boundaries  were 
located  10  chords  away  from  the  wing,  outlet  boundary  -  20  chords  away.  Side 
boundary  was  located  one  span  away  from  the  wing.  Such  choice  of  computational 
domain  is  a  compromise  between  computational  expenses  and  influence  of  boundary 
conditions.  The  computational  model  consists  from  more  than  1  200  000  cells,  the 
mesh  is  collocated  from  the  boundaries  to  the  wing. 

Components  of  the  free- stream  velocity  vector  were  imposed  at  inlet,  upper 
and  lower  boundaries,  according  to  the  values  of  airspeed  and  angle  of  attack.  All 
simulations  were  performed  for  the  airspeed  35  m/s,  angles  of  attack  were  0°,  2°,  4°, 
6°,  8°,  10°.  The  flow  was  considered  as  fully  turbulent  (Reynolds  number  Re  = 
4,6- 105),  thus  turbulence  parameters  were  also  imposed  at  these  boundaries  - 
turbulence  intensity  1%,  turbulence  length  scale  0,001  m. 
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Figure  7. 14  -  General  view  of  the  computational  domain 


Figure  7. 15  -  Computational  domain  as  it  is  seen  along  Z  axis  (in  spanwise  direction) 
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Figure  7. 16  -  Computational  domain  as  it  is  seen  along  X  axis  (forward  view) 

and  Y  axis  (view  from  above) 
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Atmospheric  pressure  of  1  •  1 05  Pa  was  specified  at  the  outlet  boundary.  To 
minimize  the  influence  of  the  outlet  boundary  it  was  considered  that  pressure 
distribution  is  not  uniform,  but  its  integral  mean  value  is  equal  to  the  prescribed 
pressure.  Zero  streamwise  gradients  were  prescribed  for  turbulence  parameters  at  the 
outlet. 

Symmetry  boundary  conditions  were  prescribed  for  side  boundaries,  i.e.  zero 
normal  flow  (along  Z-axis  for  this  case). 

No-slip  condition  was  used  for  wing  surface,  i.e.  all  components  of  the 
velocity  are  equal  to  zero. 

A  well-known  “k-e  RNG”  turbulence  model  in  combination  with  non¬ 
equilibrium  wall  function,  accounting  for  streamwise  pressure  gradient  influence  at 
the  velocity  profile  in  the  boundary  layer,  was  used  for  simulations. 

All  simulations  were  run  as  steady  state,  using  second-order-accurate 
differencing  schemes  for  all  flow  parameters. 

The  results  of  the  simulations  were  used  for  plotting  of  pressure  coefficient 
distribution  which  are  shown  at  figures  7.17-7.22. 
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Figure  7. 17  -  Pressure  coefficient  distribution  for  lower  and  upper  surface  of  the  telescope  wing, 

oc=  0° 


Figure  7. 18  -  Pressure  coefficient  distribution  for  lower  and  upper  surface  of  the  telescope  wing, 

a  =  2° 


Figure  7. 19  -  Pressure  coefficient  distribution  for  lower  and  upper  surface  of  the  telescope  wing, 

a=  4° 
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Figure  7.20  -  Pressure  coefficient  distribution  for  lower  and  upper  surface  of  the  telescope  wing, 

a=  6° 


Figure  7.21  -  Pressure  coefficient  distribution  for  lower  and  upper  surface  of  the  telescope  wing, 

a=  8° 


Figure  7.22  -  Pressure  coefficient  distribution  for  lower  and  upper  surface  of  the  telescope  wing, 

a=  10° 


Figures  7.23  and  7.24  show  flow  pattern  at  the  joint  between  two  parts  of  the 
telescope  wing  for  a  =  0°  and  a  =  10°.  It  can  be  seen  from  figure  7.24,  that  for  a  = 
10°  the  flow  is  very  complicated,  with  multiple  separations  and  recirculations. 
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Figure  7.23  -  Flow  pattern  at  the  joint  between  two  parts  of  the  telescope  wing,  a  =  0° 


Figure  7.24  -  Flow  pattern  at  the  joint  between  two  parts  of  the  telescope  wing,  a  =  10° 


By  integration  of  pressure  distribution,  spanwise  circulation  distributions, 
shown  at  figure  7.25,  were  obtained.  Local  increase  in  circulation  at  the  joint 
between  two  parts  of  the  wing  can  be  explained  by  suction  at  the  upper  wing 
surface,  caused  by  vortex  system  forming  at  the  joint.  Local  increase  of  circulation 
at  the  wing  tip  is  caused  by  suction  in  the  tip  vortex  core. 
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Figure  7.25  -  Spanwise  distribution  of  dimensionless  circulation  for 
the  telescope  wing  at  different  angles  of  attack: 

1  -  a  =  2°,  2  -  a  =  4°,  3  -  a=  6°,  4  -  a  =  8°,  5  -  a=  10° 

Aerodynamic  performance  was  also  estimated  from  the  results  of  simulations 
-  lift  coefficient  CL  vs.  angle  of  attack  (shown  at  figure  7.26)  and  drag  coefficient  CD 
vs.  angle  of  attack  (shown  at  figure  7.27).  Figure  7.28  shows  lift-to-drag  ratio  Ci/CD 
vs.  angle  of  attack,  and  figure  7.29  -  drag  polar,  i.e.  CL  vs.  CD.  It  can  be  seen  that  for 
this  range  of  angle  of  attack  the  lift  coefficient  increases  nearly  linear  with  the  angle 
of  attack.  Maximum  lift-to-drag  ratio  for  this  Reynolds  number  is  achieved  at  a  = 
8°.  Drag  due  to  lift  factor  K is  about  0,068. 
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a 


Figure  7.27  -  Drag  coefficient  vs.  angle  of  attack 
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Figure  7.28  -  Lift-to-drag  ratio  vs.  angle  of  attack 


Co 


Figure  7.29  -  Drag  polar 
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Figure  7.30  presents  comparison  of  different  calculations  for  the  telescope 
wing.  Spanwise  circulation  distribution  was  calculated  with  Star-CD  (a  =  2°)  and 
with  discrete  vortex  method,  linear  and  non-linear  (a  =  2°)  variant.  A  good  general 
agreement  can  be  seen.  The  main  difference  is  near  the  joint  between  two  parts  of 
the  wing.  Star-CD  results  show  sharp  increase  in  circulation,  while  linear  discrete 
vortex  method  give  quite  a  small  “junk”.  Non-linear  discrete  vortex  method 
overestimates  circulation  for  the  main  part  of  the  wing  comparing  with  two  other 
computational  models. 


Figure  7.30  -  Spanwise  distribution  of  dimensionless  circulation  for  the  telescope  wing 

Figures  7.31  -  7.36  present  chord  plots  of  the  pressure  coefficient  for  seven 
test  sections,  for  which  the  measurements  were  performed  in  the  wind-tunnel  tests. 
Note  that  for  the  Section  4  (located  at  the  main  part  of  the  wing  near  the  joint) 
pressure  distribution  is  sufficiently  different  from  the  other  sections.  This  is 
probably  explained  by  the  influence  of  local  separation  zones,  which  are  visible 
from  figures  7.23  and  7.24. 
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Figure  7.31  -  Simulation  results  for  a=  0 
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Figure  7.32  -  Simulation  results  for  a  =  2 
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Figure  7.33  -  Simulation  results  for  «  =  4 
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Section  1 


Section  2 


Section  5 


Section  6 


Figure  7.34  -  Simulation  results  for  a=6 
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Figure  7.35  -  Simulation  results  for  «  =  8 
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Figure  7.36  -  Simulation  results  for  a  =  10 
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8  Experimental  research 
of  telescope  wing  aerodynamics 

8.1  Description  of  test  models 
for  balance  tests  in  the  wind-tunnel 

In  order  to  analyze  aerodynamic  performance  of  the  telescope  wing  a  number 
of  test  models  were  manufactured.  Geometry  of  the  parts  of  these  models  presented 
at  figures  8.1  -  8.4. 


no 


Figure  8.1-  Geometry  of  the  nose  part  of  the  model  fuselage 


150 


CM 


Figure  8.2  -  Geometry  of  the  main  part  of  the  model  fuselage  (“clean  fuselage”  variant) 


Figure  8.3  -  Geometry  of  the  main  part  of  the  winged  model 
(“trapezoidal  wing-fuselage”  variant) 
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Figure  8.4  -  Geometry  of  the  main  part  of  the  winged  model 
(“telescope  wing-fuselage”  variant) 

Photos  of  the  model  elements  are  shown  at  figures  8.5  -  8.8. 


Figure  8.5  -  Nose  part  of  the  model  fuselage 
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Figure  8.6  -  Main  part  of  the  model  fuselage  (“clean  fuselage”  variant) 
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Figure  8.7  -  Main  part  of  the  winged  model  (“trapezoidal  wing-fuselage”  variant) 
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Figure  8.8  -  Main  part  of  the  winged  model  (“telescope  wing-fuselage”  variant) 


Experimental  program  was  comprised  of  tests  of  “clean”  fuselage  (figure  8.9), 
“trapezoidal  wing- fuselage”  combination  (figure  8.10)  and  “telescope  wing- 
fuselage”  combination  (figure  8.11). 
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Figure  8.9  -  “Clean”  fuselage  model 
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Figure  8.10  -  Model  of  “trapezoidal  wing-fuselage”  combination 
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Figure  8. 1 1  -  Model  of  “telescope  wing-fuselage”  combination 


8.2  Results  of  balance  tests 

Experimental  program  involved  four-fold  measurement  for  each  airspeed  (25, 
30,  35  and  40  m/s).  In  order  to  reduce  the  influence  of  the  model  asymmetry  all  tests 
were  performed  for  normal  and  for  inverted  (i.e.  rotated  180°  about  the  longitudinal 
axis)  position  of  the  model.  Thus  each  model  (“clear”  fuselage,  “trapezoidal  wing- 
fuselage”  and  “telescope  wing- fuselage”  combinations)  was  tested  16  times. 
Combinations  of  model  position  and  airspeed  were  picked  up  with  the  help  of 
random  number  generator. 

Experiments  were  performed  with  the  help  of  automated  system  based  on  T-3 
SSAU  low-speed  wind-tunnel.  This  wind-tunnel  is  equipped  with  inner-model  6- 
component  strain  balance.  Automated  system  sustains  constant  prescribed  airspeed 
during  the  experiment,  controls  changes  of  the  angle  of  attack,  acquires  data  from 
the  strain  balance  and  provides  real-time  low-level  data  processing. 
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Next  level  of  data  processing  includes  corrections  eliminating  the  influence  of 
different  collateral  factors. 

After  that  the  polynomial  fits  were  found  with  MATLAB  software  -  5th  order 
for  the  lift  coefficient  vs.  angle  of  attack  and  6th  order  for  the  drag  coefficient  vs. 
angle  of  attack.  At  figures  8.12  -  8.23  points  present  experimental  data  and  solid 
lines  are  polynomial  fits. 

Polynomial  fits  for  clean  fuselage  were  subtracted  from  the  polynomial  fits  for 
“trapezoidal  wing- fuselage”  and  “telescope  wing- fuselage”  combinations  in  order  to 
obtain  aerodynamic  performance  of  the  “isolated”  wings.  These  data  include  some 
effects  resulting  from  the  interference  between  wing  and  fuselage,  so  they  will  be 
different  from  performance  of  “clean”  wings.  But  for  the  considered  two  wings 
(trapezoidal  and  telescope)  such  effects  are  nearly  the  same,  so  the  aerodynamic 
performance  of  “isolated”  wing  can  be  used  for  comparison.  Results  are  shown  at 
figures  8.24  and  8.25.  It  can  be  noticed  that  the  influence  of  the  airspeed  is  rather 
small. 

Figures  8.26  and  8.27  and  present  drag  polars  of  the  “isolated”  wings  and 
table  8.1  presents  maximum  aerodynamic  efficiency  for  these  wings.  Comparison  of 
the  tables  \2  and  8.1  reveals  some  difference  between  calculated  and  experimental 
values  of  maximum  aerodynamic  efficiency.  This  can  be  caused  by  the  effects 
arising  from  low  Reynolds  numbers  and  from  interference  between  wing  and 
fuselage,  which  were  not  accounted  in  the  calculations.  But  the  experiments  also 
approved  that  aerodynamic  efficiency  of  telescope  wing  is  nearly  the  same  as  of  the 
trapezoidal  wing. 


Table  8.1  -  Experimental  values  of  maximum  aerodynamic  efficiency 


Airspeed,  m/s 

25 

30 

35 

40 

Telescope  wing 

15.3 

17.0 

17.8 

19.6 

Trapezoidal  wing 

14.7 

16.4 

17.8 

20.3 
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Figure  8. 12  -  Experimental  data  and  polynomial  fit  for  “clean”  fuselage  (V=25  m/s) 


-10  -5  0  5  10  a,,  deg. 

Figure  8. 13  -  Experimental  data  and  polynomial  fit  for  “telescope  wing- fuselage”  combination 

(V=25  m/s) 
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Figure  8. 14  -  Experimental  data  and  polynomial  fit  for  “trapezoidal  wing- fuselage”  combination 

(V=25  m/s) 


Figure  8. 15  -  Experimental  data  and  polynomial  fit  for  “clean”  fuselage  (V=30  m/s) 
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Figure  8. 16  -  Experimental  data  and  polynomial  fit  for  “telescope  wing- fuselage”  combination 

(V=30  m/s) 


Figure  8. 17  -  Experimental  data  and  polynomial  fit  for  “trapezoidal  wing- fuselage”  combination 

(V=30  m/s) 
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Figure  8. 18  -  Experimental  data  and  polynomial  fit  for  “clean”  fuselage  (V=35  m/s) 


Figure  8. 19  -  Experimental  data  and  polynomial  fit  for  “telescope  wing- fuselage”  combination 

(V=35  m/s) 
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Figure  8.20  -  Experimental  data  and  polynomial  fit  for  “trapezoidal  wing-fuselage”  combination 

(V=35  m/s) 
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Figure  8.21  -  Experimental  data  and  polynomial  fit  for  “clean”  fuselage  (V=40  m/s) 
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Figure  8.22  -  Experimental  data  and  polynomial  fit  for  “telescope  wing-fuselage”  combination 

(V=40  m/s) 


Figure  8.23  -  Experimental  data  and  polynomial  fit  for  “trapezoidal  wing-fuselage”  combination 

(V=40  m/s) 
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Figure  8.24  -  Aerodynamic  coefficients  of  “isolated”  telescope  wing 


Figure  8.25  -  Aerodynamic  coefficients  of  “isolated”  trapezoidal  wing 
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Figure  8.26  -  Drag  polars  of  “isolated”  telescope  wing 
1  -  V=40  m/s;  2-35  m/s;  3-30  m/s;  4-25  m/s 


Figure  8.27  -  Drag  polars  of  “isolated”  trapezoidal  wing 
1  -  V=40  m/s;  2-35  m/s;  3-30  m/s;  4-25  m/s 
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8.3  Description  of  the  telescope  wing  test  model 
for  study  of  the  pressure  distribution 

To  study  the  pressure  distribution  at  the  telescope  wing  a  test  model  with 
NACA  0012  airfoil  and  variable  span  of  the  movable  part  was  manufactured.  The 
model  draft  is  shown  at  figure  8.28.  Dash-and-dot  lines  mark  the  sections  were 
pressure  orifices  are  located.  Each  section  has  16  pressure  orifices.  Coordinates  of 
the  orifices  are  listed  in  table  8.2.  The  numbering  is  staggered  -  orifices  with  even 
numbers  are  at  one  side  of  the  wing  and  with  odd  numbers  are  at  the  other  side  of 
the  wing.  The  first  orifice  is  at  the  leading  edge.  General  view  of  the  model  at 
different  stages  of  manufacturing  is  shown  at  photos  at  figures  8.29  and  8.30. 


Table  8.2  -  Relative  location  of  pressure  orifices  in  the  section 


Point 

1 

2 

3 

4 

5 

6 

7 

8 

X 

0 

0.025 

0.040 

0.073 

0.113 

0.165 

0.225 

0.295 

Point 

9 

10 

11 

12 

13 

14 

15 

16 

X 

0.373 

0.450 

0.555 

0.650 

0.733 

0.805 

0.853 

0.908 

Dash-and-dot  lines  mark  sections  with  pressure  orifices. 
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Figure  8.29  -  Photo  of  the  telescope  wing  model  with  pressure  orifices 
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Figure  8.30  -  Photo  of  the  telescope  wing  model  with  pressure  orifices 


The  wing  model  was  fixed  at  the  turntable,  located  at  the  supporting  vertical 
plate  made  from  organic  glass.  This  allows  changing  the  angle  of  attack  of  the  wing 
(figure  8.31).  To  reduce  the  disturbances  the  leading  edge  of  the  supporting  vertical 
plate  was  sharpened  and  shaped  as  a  wedge.  The  flat  inner  surface  of  the  plate  serves 
as  a  screen,  so  the  experiment  can  be  performed  with  the  half-model  of  the  wing. 
This  allowed  to  increase  the  size  of  the  model  and  place  inside  if  it  102  polyethylene 
tubes  connecting  the  pressure  orifices  with  pressure  sensors  (figure  8.32).  The 
general  view  of  the  model  installed  in  the  test  section  of  the  wind-tunnel  is  shown  at 
figures  8.33  and  8.34. 

16  Honeywell™  Sursence®  pressure  sensors,  purchased  with  financing  of  the 
present  project,  were  mounted  at  one  board  (figure  8.35).  These  sensors  are 
connected  with  the  common  bus  to  the  data  acquisition  board.  PC  with  data 
acquisition  board  together  with  pressure  sensors  assembly  makes  automated 
pressure  measurement  system.  Developed  software  scans  the  pressure  sensors, 
processes  the  gathered  data  and  displays  the  data  at  the  computer  screen.  This 
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automated  pressure  measurement  system  sufficiently  decreases  the  time  required  for 
the  experiment  and  increases  the  reliability  of  the  data  by  averaging  a  large  number 
of  observations  for  each  pressure  orifice. 


Figure  8.31  -  General  view  of  the  model,  turntable  and  supporting  vertical  plate 
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Figure  8.32  -  General  view  of  the  tubes,  connecting  pressure  orifices  with  pressure  sensors 
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Figure  8.33  -  Location  of  the  model  in  the  test  section  of  the  wind-tunnel 


J 
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Figure  8.34  -  Location  of  the  model  in  the  test  section  of  the  wind-tunnel,  turntable  with  limb  and 

tubes  are  visible  in  the  right  part 
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Figure  8.35  -  Assembly  of  16  pressure  sensor.  Single  pressure  sensor  is  at  the  front 


8.4  Experimental  results  of  pressure  distribution  study 

Experiments  were  performed  with  the  airspeed  37  m/s.  Chord  plots  of 
pressure  coefficient  for  all  7  sections  for  angles  of  attack  a=5°  and  oc=10°  are  shown 
at  figures  8.36  and  8.37.  Comparison  of  these  experimental  data  with  the  results  of 
numerical  simulations  shown  at  figures  7.31  -  7.36  reveals  a  good  qualitative 
agreement.  Both  simulation  and  experiment  show  the  distinct  features  of  pressure 
distribution  in  section  4  and  7  of  the  model. 

Pressure  distributions  were  integrated  to  calculate  section  lift  coefficients. 
These  were  used  to  calculate  spanwise  load  distribution.  Figure  8.38  presents 
comparison  of  experimental  spanwise  load  distribution  for  a=5°  with  calculated 
load  distributions  shown  before  at  figure  7.30.  Experimental  results  displayed  with 
circles.  A  specific  “kink”  at  the  spanwise  load  distribution  at  the  joint  between  two 
parts  of  the  wing  is  captured  in  both  experimental  and  computational  study. 
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A  good  agreement  was  found  between  the  results  of  computational  and 
experimental  research  of  the  telescope  wing  aerodynamics.  This  can  be  considered 
as  successful  cross-validation  of  the  computational  methods  and  software  (including 
specially  developed  in  SSAU)  and  experimental  techniques  and  methodology. 

The  revealed  flow  features  at  the  joint  between  two  parts  of  the  wing  can  be  a 
subject  of  further  research  aimed  at  flow  control  in  this  area  in  order  to  increase 
aerodynamic  performance  of  the  telescope  wing. 
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Section  7 

Figure  8.36  -  Chord  plots  of  pressure  coefficient,  a  =  5° 
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Section  6 


Figure  8.37  -  Chord  plots  of  pressure  coefficient,  a  =  10° 
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Figure  8.38  -  Comparison  of  experimental  (a  =5°) 
and  calculated  spanwise  load  distributions  for  the  telescope  wing 
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9  Research  of  telescope  wing  weight 
and  aerodynamic  efficiency 

9.1  Preliminary  analysis 

For  the  tasks  of  the  present  work  a  following  algorithm  was  proposed: 

1.  Create  two  solid  bodies  -  fixed  and  mobile  parts  of  a  wing  (figure  9.1). 

2.  Assume  that  these  "solids"  can  have  nonuniform  density  distribution. 

3.  Join  these  solids  in  different  ways  -  Gevers'  version  (#1)  [7]  and  so  on 
(figure  9.2). 

4.  Optimize  the  density  distribution  within  these  solid. 

5.  Calculate  LCF  and  new  dimensionless  criteria  of  weight  efficiency. 


Figure  9.1-  Fixed  and  mobile  parts  of  a  telescope  wing 


For  optimization  of  density  distribution  in  item  5  of  the  algorithm  for  LCF 
calculation  an  iterative  procedure  is  used,  which  operates  as  following: 

1.  Some  initial  density  distribution  po/,  usually  uniform,  is  assigned  for  the 
whole  structure. 

2.  Stresses  within  all  elements  are  calculated  and  for  each  element  the 
equivalent  stress  a  e0^v  is  defined. 
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3.  New  element  densities  p/(  are  calculated  as: 


Pi;  = 


o 


eqv 


[a] 


(9.1) 


where  [<T ]  is  the  allow  stress  at  unit  density. 

The  procedure  is  repeated  until  stabilization  of  density  distribution. 

At  each  iteration  after  step  2  the  value  of  LCF  is  calculated  in  order  to  monitor 
its  change  during  the  optimization  procedure. 


© 


Rigid  connection 


Figure  9.2  -  Schemes  of  wing  parts  joints 

FEM  model  shown  at  figure  9.3  was  developed  for  weight  efficiency  analysis 
of  telescope  wing. 

Geometry  of  model  #1  of  the  telescope  wing  is  shown  at  figure  9.4. 
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Figure  9.3  -  Telescope  wing  FEM  model 
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Figure  9.4  -  Geometry  of  model  #1  of  the  telescope  wing 


To  estimate  weight  and  aerodynamic  efficiency  of  telescope  wing  the  same 
analysis  was  carried  out  for  equivalent  conventional  trapezoidal  wing  with  the  same 
area,  span,  taper  ratio  X  =  0.5  and  dimensionless  airfoil  thickness  T  =  12%. 
Geometry  of  this  wing  and  its  FEM  model  is  show  at  figures  9.6  and  9.5. 
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Figure  9.5  -  Equivalent  trapezoidal  wing  FEM  model 


Figure  9.6  -  Geometry  of  the  equivalent  trapezoidal  wing 


For  preliminary  weight  calculations  uniform  load  distribution  was  used.  The 
following  results  were  obtained: 

1.  Comparative  analysis  of  three  different  designs  of  the  joint  between  fixed 
and  movable  (extendable)  parts  of  the  telescope  wing  have  shown  that  among  them 
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the  most  preferred  is  design  #2,  presented  at  figure  9.2.  It  has  15%  better  weight 
performance  than  design  #1  (Gevers)  [7]. 

2.  For  the  uniform  wing  load  distribution  the  calculations  of  the  dimensionless 
coefficient  of  the  load-carrying  factor  yield  the  following  results: 

for  the  telescope  wing  shown  at  figure  9.4  Ck  =  31,1; 

for  conventional  equivalent  trapezoidal  wing  with  the  same  aspect  ratio  and 
area,  and  with  taper  ratio  X  =  0.5,  as  shown  at  figure  9.6  CK  =  30,86. 

Thus  an  interesting  result  was  obtained.  Telescope  wing  can  has  nearly  the 
same  weight  as  conventional  trapezoidal  wing  with  equivalent  geometry.  Such 
unexpected  result  can  be  explained  by  the  following  reasons: 

First,  the  design  #2  of  the  joint,  which  is  shown  at  figure  9.2,  does  not  cause 
sufficient  increase  of  mass. 

Second,  for  the  assumed  uniform  wing  load  distribution  the  movable 
(extendable)  part  of  the  telescope  wing  is  subjected  to  less  net  load  than  the 
corresponding  part  of  conventional  wing. 

That  is  why  in  the  following  calculations  with  the  account  for  spanwise  load 
distribution  the  design  #2  was  used  for  the  joint  (figure  9.2). 
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9.2  Weight  and  aerodynamic  analysis  of  telescope  wings  with 
the  account  for  spanwise  load  distribution 

Three  pairs  of  wings  were  considered  (figure  9.7).  An  equivalent  trapezoidal 


Figure  9.7  -  Geometric  characteristics  of  considered  wings 
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wing  with  the  same  area,  span  and  taper  ratio  X  =  0.5  corresponds  to  each  telescope 
wing.  Airfoils  with  dimensionless  thickness  T  =  12%  were  used  for  all  calculations. 

Table  9.1  presents  the  aerodynamic  performance  of  these  wings,  showing 
their  high  aerodynamic  efficiency. 


Table  9.1  -  The  aerodynamic  performance  of  considered  wings  (see  section  7) 


Wing 

0 

A 

c 

Wa 

K 

cDo 

L/Dm 

la 

0.75 

9.818 

4.9197 

0.0354 

0.00636 

33.290 

lb 

0.75 

10.058 

4.9362 

0.0323 

0.00643 

34.692 

2a 

0.50 

10.800 

5.0183 

0.0301 

0.00639 

36.051 

2b 

0.50 

10.941 

5.0238 

0.0297 

0.00653 

35.895 

3a 

0.25 

12.000 

4.9096 

0.0311 

0.00638 

35.461 

3b 

0.25 

12.000 

5.1146 

0.0272 

0.00663 

37.257 

Figure  9.20  explains  the  value  of  ct . 

For  weight  analysis,  which  is  discussed  further,  spanwise  load  distributions 
shown  at  figures  9.8  -  9.1 1  were  used. 


Figure  9.8  -  Wings  la,  lb  dimensionless  velocity  circulation 
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0  0.2  0.4  0.6  0.8  z 

Figure  9.9  -  Wings  2a,  2b  dimensionless  velocity  circulation 


Figure  9.10  -  Wings  3a,  3b  dimensionless  velocity  circulation 
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0  0,2  0,4  0,6  0,8  z 


Figure  9. 1 1  -  Wings  la  -  3a  dimensionless  velocity  circulation 
FEM  models  of  the  wings  2a  and  2b  with  loads  are  shown  at  figures  9.12  and 

9.13. 

Figures  9.14  and  9.15  illustrate  the  proposed  load-carrying  scheme  and  load 
transfer  in  the  joint  of  telescope  wing. 

Results  of  weight  analysis  for  all  of  the  six  wings  are  presented  in  table  9.2. 
Columns  depict  the  change  of  CK  during  the  process  of  material  distribution 
optimization. 


Figure  9. 12  -  Loads  via  dimensionless  velocity  circulation  for  a  telescope  wing 
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z 


Figure  9. 13  -  Loads  via  dimensionless  velocity  circulation  for  equivalent  trapezoidal  wing 

Calculations  of  CK  for  wings  3a  and  3b  gave  the  values  23.39  for  telescope 
wing  and  28.98  for  equivalent  trapezoidal  wing.  This  yields  that  telescope  wing  can 
theoretically  have  even  less  weight  than  aerodynamically  similar  trapezoidal  wing. 
Prima  facie,  this  result  seems  paradoxical  and  requires  explanation. 


Figure  9. 14  -  Scheme  of  contact  elements  displacements 
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Figure  9. 15  -  Contact  elements  loading 


The  weight  of  high-aspect -ratio  wing  is  mostly  influenced  by  the  magnitude 
of  bending  moment  and  airfoil  thikness  of  the  wing.  Figure  9.16  shows  the 
distribution  of  bending  moment  Mx  for  wings  3a  and  3b.  It  can  be  seen  that  the 
telescope  wing  has  considerably  smaller  bending  moment  along  its  whole  span.  The 
weight  of  unit  length  of  the  high-aspect -ratio  wing  can  be  considered  proportional  to 
M x/trmx  in  the  given  cross-section.  Thus  it  can  be  taken  as  the  first  estimation  that 
the  weights  of  the  considered  wings  is  determined  by  two  controversial  factors  Mx 
and  tmax.  Distributions  of  Mxjtmsx  for  wings  3a  and  3b  are  shown  at  figure  9.17. 
This  plot  explains  why  the  telescope  wing  can  be  lighter  than  trapezoidal. 
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Table  9.2  -  Calculation  of  Ck  criteria  for  wings 


Name 

Wing  la 

Wing  1b 

Wing  2a 

Wing  2b 

Wing  3a 

Wing  3b 

Cki 

31,9970 

31,5810 

33,5743 

35,7660 

32,8615 

40,9904 

Ck2 

29,4924 

27,3511 

30,2403 

30,9658 

28,8096 

35,4783 

Ck3 

28,3572 

25,4112 

28,5571 

28,7560 

26,8196 

32,9297 

Ck4 

27,6988 

24,2965 

27,5169 

27,4818 

25,6503 

31 ,4570 

Ck5 

27,2615 

23,6102 

26,8121 

26,6946 

24,9175 

30,5437 

Ck6 

26,9466 

23,1572 

26,3070 

26,1767 

24,4245 

29,9438 

Ck7 

- 

22,8375 

25,9323 

25,8114 

24,0761 

29,5208 

Ck8 

- 

22,6049 

25,6476 

25,5439 

23,8219 

29,2110 

Ck9 

- 

22,4331 

- 

25,3466 

23,6343 

28,9811 

Ckio 

- 

- 

- 

- 

23,4933 

- 

Ckii 

- 

- 

- 

- 

23,3864 

- 

8,50E+07 

8,00E+07 

7,50E+07 

7,00E+07 

6,50E+07 

6,00E+07 

5,50E+07 

5,00E+07 

4,50E+07 

4,00E+07 

3,50E+07 

3,00E+07 

2,50E+07 

2,00E+07 

1.50E+07 

1.00E+07 

5,00E+06 

0,00E+00 


Figure  9. 16  -  Comparison  Mx  for  wings  3a  and  3b 
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Figure  9. 17  -  Comparison  Mx/tmax  for  wings  3  a  and  3b 

According  to  the  results  of  our  research  for  the  telescope  wing  the  joint  itself 
doesn’t  yield  sufficient  increase  of  CK.  We  were  trying  to  intentionally  worsen  and 
place  contact  pairs  in  the  area  of  ?max  of  the  main  part  of  the  wing.  This  leads  to 
greater  load  on  the  tip  rib.  Sufficient  increase  of  CK  by  2.3%  was  obtained  only  for 
the  wing  la.  It  can  be  assumed  that  the  presence  of  the  joint  causes  only  some 
increase  in  constructional  coefficient  (p  (full  mass  coefficient). 


9.3  Comparison  weight  and  aerodynamic  efficiency 
of  the  telescope  and  eqvivalent  trapezoidal  wings 

Research  results  for  weight  and  maximum  lift-to-drag  ratio  L/Dmax  of  telescope 
and  equivalent  trapezoidal  wings  are  plotted  at  figures  9.18,  9.19. 

Figure  9.20  presents  all  the  necessary  geometrical  explanations. 

From  figures  9.18  and  9.19  it  follows  that  the  considered  telescope  and 
equivalent  trapezoidal  wings  have  nearly  the  same  weight  and  aerodynamic 
efficiency.  (Dimensionless  coefficient  CK  determines  the  wing  weight  with  an 
accuracy  to  constant  multiplier.  See  section  4.4,  (4.17,  4.18,  4.19).) 

For  some  geometrical  parameters  the  telescope  wing  theoretically  can  have  less 
weight  then  the  equivalent  trapezoidal  wing. 
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Figure  9. 18  -  Comparison  Ck  for  wings  after  iterations 
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Figure  9. 19  -  Comparison  maximum  aerodinamic  efficiency  for  wings 

Addition .  Sign  "o"  on  figure  9.18  presents  a  load-carrying  factor  coefficient 
calculated  for  uniform  load  distribution.  This  results  shows  the  simplified  load 
application  gives  ability  to  adequate  comparison  of  various  configurations  of  wings. 
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Figure  9.20  -  Explanation  for  choise  of  geometric  characteristics  equvalent  trapezoidal  wing 
c  . 

Ct  = - dimensionless  chord  of  rectangular  outer  wing  part 

C  root 

(referenced  to  the  tip  chord  of  the  main  part  of  the  wing); 
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10  Results  and  conclusions 

1.  Accuracy  analysis  for  some  known  wing  weight  formulas  was  performed, 
showing  that  they  give  too  wide  spreading  of  the  results  (tens  of  percents)  and  are 
not  able  to  account  for  some  features  of  morphing  wing. 

2.  New  method  was  proposed  in  order  to  increase  accuracy  of  weight  design  and 
provide  adequate  comparison  of  different  engineering  solutions.  In  principal  new  in 
this  method  are  use  3D -models  with  optimization  of  variable  density  of  elements  for 
calculation  of  load-carrying  factor  for  theoretically  optimal  structures,  and 
introduction  of  dimensionless  load-carrying  factor  coefficient,  showing  the  weight 
perfection  of  load-carrying  structure.  This  criteria  allows  objective  comparison  of 
different  structure,  thus  creating  knowledge-  and  database  of  different  engineering 
solutions.  New  weight  formulas  developed  basing  on  this  criteria,  which  allow  to 
calculate  absolute  wing  mass  and  wing  mass  fracture.  The  proposed  approach  do  not 
impose  any  limitations  for  external  shapes  and  boundary  conditions  of  wing 
structures. 

3.  For  telescope  wings  a  wide  range  of  aerodynamic  research  was  performed. 
Spanwise  load  distribution  was  calculated  with  three  different  computational 
methods.  It  was  shown  that  discrete  vortex  method  provides  reasonably  accurate 
results.  This  method  was  used  for  calculation  of  spanwise  load  distribution  for 
telescope  and  trapezoidal  wings. 

4.  Three  test  models  for  wind-tunnel  testing  were  designed  and  built  in  order  to 
verify  computational  methods  used  for  telescope  wing  aerodynamic  calculations. 
Two  models  for  load  tests  and  one,  with  pressure  orifices,  for  study  of  pressure 
distribution  in  the  joint  area  of  the  telescope  wing.  Wind-tunnel  testing  verified  the 
possibility  of  obtaining  high  aerodynamic  efficiency  of  telescope  wings  and  revealed 
some  special  features  in  the  flow  pattern  in  the  joint  area. 

5.  The  rational  type  of  the  telescope  wing  joint  was  determined  as  a  result  of 
weight  analysis.  The  research  of  telescope  wings  have  shown  that  they  can  have 
nearly  the  same  weight  efficiency  and  aerodynamic  performance  as  equivalent 
trapezoidal  wings  of  the  same  span  and  area. 

6.  Preliminary  results  of  telescope  wing  optimisation  have  shown  that  high 
weight  efficiency  and  aerodynamic  performance  are  obtained  with  small  relative 
(referency)  chord  of  the  movable  part  of  the  wing,  about  0.25.  Such  variant  of  the 
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telescope  wing  can  be  considered  as  perspective  for  liner  aircraft  and  other 
applications. 

7.  Comparison  of  the  results  of  weight  estimations  for  telescope  wings,  obtained 
with  two  different  approaches  -  with  uniform  load  distribution  and  with  spanwise 
load  distribution  calculated  with  discrete  vortex  method  -  has  shown  that  the  first 
approach  can  be  used  for  comparing  weight  efficiency  of  different  designs.  But  for 
weight  design  the  second  approach  should  be  used,  because  the  difference  can  make 
up  to  10-15%. 

8.  Calculations  of  full  mass  coefficient  for  maneuvarable  aircraft  wing  with  large 
area  of  movable  parts,  and  for  flap  and  slat  of  liner  aircraft  resulted  in  a  wide 
spreading  of  values,  which  points  to  the  necessity  of  additional  research  of  this 
coefficient,  concerning  mainly  the  account  for  multiple  load  cases. 

9.  Comparison  of  numerical  simulations  of  flow  past  the  model  telescope  wing 
with  the  experimental  results  from  pressure  orifices  model  have  shown  a  good 
agreement.  Thus  a  mutual  verification  of  computational  and  experimental 
methodology  was  performed.  The  technology  used  for  creation  of  pressure  orifices 
model  and  pressure  distribution  measurement  with  Honeywell™  Sursence® 
DUXL10D  pressure  sensors  with  automated  data  processing  has  proved  to  be  highly 
reliable  and  accurate  and  can  be  used  further  for  quick  research  of  the  same 
problems. 

10.  The  present  research  have  required  coordinative  work  of  the  specialists  in 
different  fields  -  computational  modelling  and  optimization  of  the  structures, 
computational  and  experimental  aerodynamics.  For  the  purpose  of  this  research 
necessary  licensed  software  and  experimental  equipment  was  purchased,  and  some 
additional  software  was  developed.  The  important  result  of  this  project  is  the 
creation  of  the  team  which  is  capable  for  solving  wide  range  of  problems  dealing 
with  morphing  wings  and  similar  research. 
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11  Future  work  recommended 


1 .  Development  of  the  proposed  method  for  weight  analysis  in  order  to  account 
for  multiple  load  cases. 

2.  Further  multi-disciplinary  research  of  morphing  wing  for  the  purposes  of 
optimization  of  their  parameters  and  search  of  their  efficient  applications. 


REO-1386  Research  of  Morphing  Wing  Efficiency 

Final  Technical  Report  /  Samara  State  Aerospace  University,  AVIKON  Institute 
1  October  2003  -  1  June  2004 


Ill 


12  Plan  of  publishing 


Preparation  of  the  following  papers  is  proposed: 

1.  Application  of  3D  modelling  to  weight  estimations  of  aircraft  structures. 

2.  Features  of  telescope  wing  aerodynamics. 

3.  Estimation  of  weight  and  aerodynamic  performance  of  telescope  wings. 
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